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Abstract 

The time-dependent thermopower is analyzed through an interacting 
quantum dot coupled to a time-dependent gate voltage and under the 
influence of an external magnetic field using the Keldysh nonequilibrium 
Green's function formalism. Formal expressions of the electrical and ther- 
mal conductances, thermopower, and thermoelectrical figure of merit are 
obtained. The influence of the magnetic field on the displacement current 
and the heat current is studied. Results show that although applying 
time-dependent gate voltage results in the enhancement of the Seebeck 
coefficient, the electron-electron interaction gives rise to a significant re- 
duction in the thermopower. The reason for why applying time dependent 
gate voltage results in the enhancement of the thermopower is also ana- 
lyzed. 

1 Introduction 

Thermopower - the ratio of induced voltage to an applied temperature gradient 
across the sample at the state of vanishing current - is one of the oldest issues of 
solid-state physics. The study of thermopower in nanoscale devices has attracted 
a lot of attention in recent years because of recent developments in fabricating 
and utilizing them. In addition, deviations from the Wiedmann-Franz law in 
nanostructures result in significant enhancements in the thermopower of such 
devices [JJ . Various models have been suggested for the study of thermoelectric 
transport through quantum dots (QDs) and molecular junctions [U [3J HJ [51 EH H] i 
nanotubes and nanowires [HI HJ [TU] , strongly correlated nanostructures [TTJ [T^l 
I13j . etc. In addition, measurement of the thermopower through nanostructures 
has been an interesting topic in recent years [14j [T5l [TU [TTJ [18] . 

Recently, Crepieux and co-workers [19] proposed applying time-dependent 
gate voltage results in the enhancement of the thermopower. Their results 
showed that the thermopower is enhanced by up to 40 %. Time-dependent 
transport has been extensively studied both theoretically and experimentally 
[2ni[nj[21[23[Ml[2Sl[2Sl[2I]- However, the time-dependent heat current and 
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the thermopower are novel phenomena requiring more attention. In this paper, 
we consider an interacting QD coupled to metal leads. Although the system 
seems to be similar to what was considered in Ref. [T9] , the effects of electron- 
electron interactions and the magnetic field are also taken into account. Using 
the Keldysh nonequilibrium Green's function formalism [28, 25], expressions for 
the electrical and thermal conductance are obtained. The results show that 
the time-dependent thermopower is significantly reduced by electron-electron 
interactions. The influence of the magnetic field on the heat current is also 
investigated. 

In the next section, the heat current is evaluated using the nonequilibrium 
Green's function formalism. We use the Hartrce approximation and hence the 
heat current is related to the averaged electron density. In sect. 3, numerical 
results arc presented and in the end, conclusion is presented. 



2 Model 

We consider an interacting QD coupled to metal leads and under the influence 
of a step-like gate voltage pulse. The Hamiltonian describing the system is given 
as follows: 

H = Y1 £ »k*cl kr7 c aka + £a{t)n a + U n t ni + ^[V^c^c^ + H.C] (1) 

aka <y aka 

where c a ka(c ak(7 ) destroys (creates) an electron with wave vector k, spin a, and 
energy £ a ka in lead q ( q = L or R ). is the annihilation (creation) 

operator for the dot and n a = S a d a is the occupation number. £ a (t) denotes 
the time-dependent energy level of the QD defined as e a (t) = e" + A d Q(t) 
where is the time-variation of the gate voltage. e° = £o ± E z (plus sign 
for spin-up) is the time-independent energy level of the QD and E z is Zeeman 
splitting induced by an external magnetic field. U and V a ka stand for Coulomb 
repulsion and tunneling strength between the dot and the lead a, respectively. 
The time-dependent heat current (J h ) is obtained from the difference between 
the energy current (I e ) and the charge current (I q ) according to 



/ h (t)=7 e (<) _M!) 



TO (2) 



where 1% — J2kaea < dt C aka Ca ka- >, la ~ ^2kaea £ak <? < dt C aka Cak a >, 

and fj, a denotes the chemical potential of lead a. The Keldysh nonequilibrium 
Green's function formalism is used to obtain the energy current. It is straight- 
forward to show that by means of the Green's function of the isolated leads, the 
energy current is given as {h = 1) [19] 

I c a (t) - 2i?e |V QfcCT | 2 / dtikatoe- 16 -*'^-*) (3) 

[Ue aka )Gl a (tM) + G< a {tM)] 
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where / a (e) = [1 + exp((e — fj, a ) / 'fcT a )] _1 is the Fermi distribution function 
of the a th lead and T a denotes the temperature of the lead. G r aa ,{t,t') and 
G'S , (i, t') are the retarded and the lesser Green's functions of the interacting 
QD, respectively. 

The Green's function of the QD is obtained from a Dayson-like equation as 



Gw(t,t') = gw/(r, r') + / dt 1 dt 2 g a a' (r, ri)S CT (n, T 2 )G crcr ' (t 2 , r) 



(4) 



where C denotes the contour integral, and E CT (ii,i2) = X^afc IVafco-Pffafco-Cii, £2) 
is the self-energy operator in which g a ka(t,t') is the isolated lead Green's func- 
tion. g aa i is the isolated QD Green's function satisfying [3D] 

[i^ -e a {t)-U< n s {t) >}g aa ,{t,t') = 5{t - i')<W (5) 

For deriving the above equation, the decoupling approximation introduced in 
Ref. [21] has been used, i.e., < {n s (t)d a (t),dl,(t')} >«< n s (t) > g acr i (t,t'). < 
n a (t) > denotes the time averaged electron density with spin er, and a is opposite 
of a. This approximation is reasonable under conditions that the temperature 
is low enough (lower than the level spacing) and the bias voltage is small. In the 
following, we change the sum over k into an energy integral and use the wide 
band approximation, i.e., J2k \ Vaka\ 2 = J defl-xT^ where = 2Trp a \V a k a \ 2 
is the spin-dependent tunneling rate. The retarded and lesser Green functions 
are obtained from Eq. 4 by means of the Langreth continuation theorem |31j . 
Then, the time-dependent heat current is given as 

^) = -E r «E r «'/ ^-^)/a'(£)K<(M)| 2 (6) 

<J a' 

+ J ^(e-^)f a (e)hn{AZ(e,t)}] 
where A°(e,t) is given by [20] 

[ £ - e na + i/2r CT ] - ^[e~e^-^ d +i/2T-]t 

a l£ ' ' [e- e na + i/2T*] [e - e na - A d + i/2T°] 1 ' 

where e na = e a + U < ns > and T a = T a L + T^,. It is straightforward to show 
the time averaged electron density is obtained from [20] 

<n CT (t)> = T lim o ^^r- J ^-Ue) jf dt\A° a {e, t)| 2 (8) 



E 



[(e - e na f + (if )2][( £ - £na A d f + 



Equation 8 should be solved self-consistently to obtain the electron density. In 
the following, we assume that = r and use T as the energy unit [32] and 
h/To as the time unit. 
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Figure 1: (a) 1% versus time for = 0.5 (solid), = 1 (dashed) and A^ = 
1.5 (dashed-dotted). (b) and (c) show spin- resolved displacement current for 
E z = (solid), E z — 1 (dashed) and E z = 2 (dashed-dotted). Parameters are 
Ti = 2,T 2 = 0.1, £ = 0,U = 2, and n L(R) = +(-)2.5. For (b) and (c), we set 
A d = l. 
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3 Results and Discussion 



Figure 1 shows the spin-dependent displacement current (I%(t) — I£(t) + 
as a function of the amplitude of the gate voltage (A<j) (Fig. la) or the magnetic 
field (Figs, lb and c). The displacement current describes the time evolution 
of the electron density of the QD. As we expect, /a is zero at t < 0, because 
the system is in the steady state. Upon applying the gate voltage , a signifi- 
cant reduction in I£ is observed when < t < 2. The gate voltage shifts the 
QD's energy level toward the chemical potential of the emitter (left lead) and 
hence, a sudden decrease in the population of the charge occurs. The reduction 
of the population leads to J£ < 0. During the time, the system approaches 
the new steady state and, therefore, 1% becomes zero again. This behavior of 
the displacement current was recently reported in Ref. [37] • Unlike the results 
obtained in Ref. [33], the displacement current does not show any fluctuation. 
Indeed, in Ref. [33] the coupling at t = adds some energy to the system but 
here, the gate voltage just changes the position of energy levels. Figures lb and 
c describe the behavior of the spin-up and spin-down displacement currents in 
the response to the external magnetic field, respectively. It is observed that /I 
is not sensitive to the magnetic field. Indeed, in the presence of the magnetic 
field, the population of the spin-up is significantly reduced, because the entrance 
of the spin-up electron into the QD requires the energy to be more than \x\, due 
to electron-electron interactions. The behavior of is more interesting in the 
response to the magnetic field. Unlike the spin-up level, the spin-down level is 
always inside the bias window and the increasing magnetic field gives rise to the 
enhancement of the spin-down electron population. It is found that becomes 
positive in a strong magnetic field (E z = 2) when 1.5 < t < 2.5. It comes from 
the inequality in the temperature of the leads. 

The spin resolved heat current is plotted in Fig. 2 for different magnetic 
fields. The typical coherent oscillations in the heat current are observed. It is 
interesting to note that the frequency of the oscillations is spin- and Coulomb 
repulsion-dependent and given by ftu a — \Ef — e a — U < n s > — Ad|. The 
existence of such beats in the charge current was previously reported (22] [27] ; 
however, for the heat current more experimental results are needed. At t < 0, 
or t — > co, the heat current approaches constant values. It is also observed that 
the variations of the heat current in the right lead are more significant than 
in the left one because of Tr < 2l- Notice, the heat current becomes positive 
when the energy level of the QD is outside the bias window. Therefore, J^t > 
for E z = 1 and 2. In the case of 1^,, the magnitude of the heat current is 
enhanced by an increase of the magnetic field because of //l » £[.. The story 
is completely different about I^ a , because an increase of magnetic field leads to 
Ea- > MR an d hence, the magnitude of the heat current is decreased. 

In the linear response limit, the current is given by I(t) = Gy{t)AV + 
Grit) AT where Gyif) and Gr(t) are the electrical conductance and the thermal 
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Figure 2: Spin-dependent heat current against time. E z = (solid), E z = 1 
(dashed) and E z = 2 (dashed-dotted). Other parameters are the same as Fig. 
1. 



coefficient, respectively. Setting I{t) = 0, the Seebeck coefficient is defined as 

From Eq. 6, we obtain the time-dependent conductance coefficients as fol- 
lows [31]: 

<M*) = "E^/ ^f(e)Im(A^e,t)) (10a) 

G T {t) = ^E^P / ^-(e-e F )f(e)Im(A°(e,t)) (10b) 

where /'(e) = d/de/(e). Note that the above equations are only valid under 
conditions that AV^ and AT are small in comparison to £o and A^ — > 0. In 
these conditions, we have \A a (e,t)\ 2 = -2/r°'Im[A cr (e, t)] [19 . It is interest- 
ing to note that the imaginary part of A(e, t) is thought of as the generalized 
spectral function. The time-dependent Seebeck coefficient is plotted in Fig. 3a. 
One observes that there is a sudden increase in thermopower once the gate 
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Figure 3: a) Thermopower, (b) conductance, and (c) thermal coefficient versus 
time. Parameters are T x = T 2 = 0.1, e = 0.5, E z = 0, U = 2 and A d = 0.05 
(solid) and Ad = 0.1 (dashed). Inset shows the generalized spectral function as 
a function of energy at t = 1 (solid line) and t = 7 (dashed line). Normalized 
/'(e) is also plotted in gray. 



voltage is applied. Such a time-dependent enhancement has been recently re- 
ported in a non-interacting QD |19j . Indeed, the more sensitivity of the system 
to the temperature in the transient regime results in the enhancement of the 
thermopower. Since Ad is small, the imaginary part of A"(e,t) is composed of 
a dominant Lorentzian-like part where s na is its center, and a corrective part 
which is on the order of Ad and decays during the time according toe" 1 / 2 ^ 4 . It 
is interesting to note that the correction part plays the main role in the enhance- 
ment of the thermopower. /'(e) has a Lorentzian shape which is centered near 
the chemical potential of the lead, see inset of Fig. 3. Under these conditions, 
the corrective term causes Im(^4 CT (e, f)) to become larger at initial times after 
applying gate voltage. It leads to the enhancement of the thermopower. This 
term decays with time and, as a result, the thermopower becomes constant 
again. Indeed, the increasing generalized spectral function near the chemical 
potential of the leads results in the enhancement of the thermopower upon ap- 
plying gate voltage. The generalized spectral function and /'(e) are plotted in 
the inset of Fig. 3a. It is evident that the increase of Ad leads to more increase 
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of the generalized spectral function and, as a result, the thermopower is more 
enhanced. 

Figs. 3b and c show the time evolution of Gy and Gt as a function of Ad, 
respectively. It is observed that an increase in Ad leads to an increase of the 
variation amplitude of Gy and Gt- 

The influence of the Coulomb interaction on the thermopower is analyzed in 
Fig. 4 as a percentage of (5' max — S sat )/S sat where S" nax is the maximum value 
of the thermopower and S sat = S(t — > oo). It is found that the thermopower is 
reduced by up to 25 % for strong electron-electron interactions. This reduction 
is more significant at high A^. Therefore, the predicted enhancement of the 
thermopower up to 40 % cannot be observed in strong electron-electron interac- 
tions. Indeed, an increase of correlation between electrons results in a decrease 
of the thermopower. 




Figure 4: Percentage of (5" nax - S(t ->• oo))/S(t ->• oo). We set e = 0.2. Other 
parameters are the same as in Fig. 3. 



4 Conclusion 

In this paper, we analyze the time-dependent Seebeck coefficient through an 
interacting quantum dot subject to a magnetic field. The formal expression of 
the thermopower is obtained using the nonequilibrium Green's function formal- 
ism. The influence of the magnetic field on the displacement and heat currents 
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is examined, and it is observed that the current of the left and right leads are 
different from each other in the response to the magnetic field. Spin-dependent 
beats in the heat current are also observed. We find that the thermopower is 
reduced by up to 25 % under certain conditions. The behavior of the electrical 
and thermal conductances in the response to the time-dependent gate voltage 
are also examined. 
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